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^ Abstract. The orthogonality properties of certain subspaces as- 
sociated with bivariate Bernstein-Szego measures are considered. 
It is shown that these spaces satisfy more orthogonahty relations 
than expected from the relations that define them. The results are 
CSl used to prove a Christoffel-Darboux like formula for these mea- 
sures. 

> 

u 

1. Introduction 

^ In the study of bivariate polynomials orthogonal on the bi-circle 

progress has recently been made in understanding these polynomials 
T— I in the case when the orthogonality measure is purely absolutely con- 

^ tinuous with respect to Lebesgue measure of the form 

where Pn,miz,w) is of degree n in z and m in w and is stable i.e. is 
^ nonzero for \z\, \ w\ < 1 and da is the normalized Lebesgue measure on 

the torus T^. Such measures have come to be called Bernstein-Szego 
measures and they played an important role in the extension of the 
^ Fejer-Riesz factorization lemma to two variables [1], [2], [5]. In 

^ particular in order to determine whether a positive trigonometric poly- 

nomial can be factored as a magnitude square of a stable polynomial 
an important role was played by a bivariate analog of the Christoffel- 
Darboux formula. The derivation of this formula was non trivial even 
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if one begins with the stable polynomial Pn,m, P, [3], [1], [9]. This for- 
mula was shown to be a special case of the formula derived by Cole and 
Wermer [1] through operator theoretic methods. Here we give an al- 
ternative derivation of the Christoffel-Darboux formula beginning with 
the stable polynomial Pn,m- This is accomplished by examining the 
orthogonality properties of the polynomial Pn,m in the space L'^{dfi). 
These orthogonality properties imbue certain subspaces of L?{d^) with 
many more orthogonality relations than would appear by just examin- 
ing the defining relations for these spaces. 

We proceed as follows. In section 2 we introduce the notation to be 
used throughout the paper and examine the orthogonality properties 
of the stable polynomial p„,m in the space L'^{dfi). We also list the 
properties of a sequence of polynomials closely associated with Pn,m- 
In section 3 we state, and in section 4, prove, one of the main results 
of the paper on the orthogonality of certain subspaces of L'^{dfi). We 
also establish several follow-up results which are then used in section 
5 to derive the Christoffel-Darboux formula. The proof is reminiscent 
of that given in [3] and f6]. In section 6, we study connections to the 
parametric moment problem. 

2. Preliminaries 

Let pn,m £ C[z, w] be stable with degree n in z and m in w. We will 
frequently use the following partial order on pairs of integers: 

{k,l) < iS k < i and / < j. 

The notations ^, ^ refer to the negations of the above partial order. 
Define 

Pn,miz,w) = Z''w"'pn,m{^/Z,l/W). 

When we refer to "orthogonalities," we shall always mean orthogo- 
nalities in the inner product (■, ■) of the Hilbert space L'^{l/\pn,m\'^da) 
on T^. Notice that L'^{l/\pn,m\'^da) is topologically isomorphic to L^(T^) 
but we use the different geometry to study Pn,m- 

The polynomial Pm,n is orthogonal to more monomials than the one 
variable theory might initially suggest. More precisely. 

Lemma 2.1. In L^{l/\pn,m\'^da), pn,m is orthogonal to the set 

{zV:(z,j)^(0,0)} 

and Pn,m is orthogonal to the set 

{z'w^ : ^ (n,m)}. 
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Figure 1. Orthogonalities of p^^^ 
j 




Pn,m -L z^w^ for in the above region 



Proof. Observe that since l/pn,m is holomorphic in W 

{z'w\pn,m)^ f Z'w^pn,m{z,w)- ^ rr^ 

7t2 \Pn,m{z,W)\^ 

= / , / a^Oif (i,j)^(0,0) 

by the mean value property (either integrating first with respect to z 
or w depending on whether i > or j > 0). The claim about Pn,m 
follows from the observation {z'^w^ ,pn,m) = {Pn,m, z"'~'^w"^'~^) . □ 

Write Pn,m{Z: W) = J2ZoPi(^)'^'- 
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Since Pn,m{z,w) is stable it follows from the Schur-Cohn test for 
stability [1] that the m x m matrix 

Ypo{l/z) Pi{l/z) 



Po[z) 



J)m-l[Z) ■■■ 
Pm{l/z) 

Ipiil/z) 



Po[z) 



L O 

o ■ 

pm{l/z)_ 



Pm[Z) 

. o 



•• Po{l/z) . 



Pm{z). 



{21] 



is positive definite for |z| = 1. Here Pj{z) = Pj{z). 

Define the following parametrized version of a one variable Christoffel- 
Darboux kernel 



Tl \_ nPn,miz,w)pn,mil/z,r]) {Z,w)pn,mil/Z,T]) 

J^[z,w,ri) — z i^z.zj 



1 — wr] 



z-[l,...,w"^-']TUz)[l,. 



m— l]t 



m—1 

= ^aj{z,w)f]\ 

j=0 

where aj{z, w), j = 0, . . . , m—1 are polynomials in {z, w), as the follow- 
ing lemma shows in addition to several other important observations. 

Lemma 2.2. Let Pn,m{z,w) be a stable polynomial of degree {n,m). 
Then, 

(1) L is a polynomial of degree {2n,m — 1) in {z,w) and a polyno- 
mial of degree m — 1 in fj. 

(2) ■; rj) spans a subspace of dimension m as rj varies over C. 

(3) L is symmetric in the sense that 

L{z,w;v) = z^^{wr]r-'L{l/z,l/w-l/r]), 

so ak = flm-fc-l- 

(4) L can be written as 

L{Z, W; 7]) = Pn,m{z, w)A{z, W, T]) + Pn,m{z, w)B{z, W, T]) 

where A, B are polynomials of degree (n, m—1, m—1) in {z,w, ff). 

Proof. The numerator of L vanishes when w = 1/f/, so the factor (1 — 
wfj) divides the numerator. This gives (1). 

For (2), when \z\ = 1 use equation (2.2). Since T^iz) > for \z\ = 1, 
L{z, w; rj) spans a set of polynomials of dimension m. 

For (3), this is just a computation. 
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m 



n 



L ± z^wj for in shaded region 



Figure 2. Orthogonalities of L. See Theorem 3.1 



For (4), observe that (suppressing the dependence of p on n and m) 



^ p{z,w)p{l/z,ri) -p{z,w)p{l/z,ri) 
1 — Wfj 

/ fj"^p(z,l/f]) — 'i]"^p(z,w)\ ^. . f fi^piz.w) — fj'^piz.llff) 
^p[z,w) \ +p[z,w^' 



1 — wrj 

V 

A(z,w,ri) 



1 — wrj 

B{z,w\rf) 



(2.3) 

□ 



3. Orthogonality relations in L^(l/|p„^mp(io-) 

Our main goal is to prove that L and oq, . . . , ctm-i possess a great 
many orthogonality relations in L'^{l/\pn,rn\^da). The orthogonality 
relations of L are depicted in Figure |2] 
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Theorem 3.1. In LF'{l/\pn,rn\^d<y), each is orthogonal to the set 
Ok ={z'w^ ■.i>n,j<0} 

U {z'w^ ■ < j < mj k} 
U {z^w-' '■ i < n,j > m} 
U {z'w'' -.i^n}. 
In L^(l/|p„ mpc^cr); L{-,-;ri) is orthogonal to the set 
O ={z''w^ : i > n, j < 0} 

U {^V : i n,0 < J < m} (3.1) 
U {z'tu^ '■ i < n,j > m}. 

Note that 

Ok = {z^'w^ : < J < m, J ^k}UO, 

m—l 

o= f]Ok. 

k=0 

Corollary 3.2. In L^{l/\pn,m\'^da), the polynomial ak is uniquely de- 
termined (up to unimodular multiples) by the conditions: 

at G span{2;W : (0,0) < (i, j) < {2n,m - 1)}, 

ak : (0, 0) < (z, j) < {2n, m - 1), j ^ k} 

U {z'w'' : < i < 2n,z 7^ n}, 

and 

\W\?= [TkAe-re^')^. 



(The last fact follows from Proposition 6.1, which is not currently 
essential.) 

Remark 3.3. We emphasize that (1) each is explicitly given from 
coefficients of Pn,rm (2) each is determined by the orthogonality 
relations in Corollary 3.2 (depicted in Figure [3]), and (3) each satisfies 



the additional orthogonality relations from Theorem |3.1 One useful 
consequence of this is that the set 

{z^ak{z, w) : j E Z,0 < k < m} 

is dual to the monomials 

: j eZ,0 <k <m} 
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Fourier support of at 



m. 

-kr 



2n 



flfc _L z'w^ for (i, j) shaded 



2n 



Figure 3. The are uniquely determined by the above 
properties. See Corollary |3.2 



within in the subspace 

S = spaTi{z^w'' : j eZ,0 < k < m}. 

Namely, 

unless ji = j2 and fci = ^2- 
In particular, if f E S, then 

/ ± z^ttk implies /(j + n,k) = 0. 



(3.2) 



4. The proof of Theorem 13.11 
We begin by writing 



m— 1 m— 1 

A{z, w;ri) = ^ Aj{z, w)f]^ B{z, w;ri) = ^ Bj{z, w)f]^ . 

j=0 j=0 



Recall equation (2.2) and Lemma 2.2 item (4). By examining coef- 
ficients of f]^ in L 

Also, Aj and Bj have at most degree j in w. To see this, recall equation 
(2.3) and observe that 

A[z,w]r]) = l^PjW^——— 
j 

which shows that Aj{z,w) has degree at most j in w (i.e. powers of w 
only occur next to greater powers of f]). The same holds for B. 
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Proof of Theorem 3J_. By Lemma 2.1, Pn,m is orthogonal to 

{^V:(z,j)^(0,0)} 
and since Ak has degree at most n in z and in if, 

Pn,mAk is orthogonal to {z'w^ : ^ {n,k)}. 

Also, 

Pn,mBk is orthogonal to {z^w^ : ^ (n,m)} 



since the orthogonality relation for Pn,m (also from Lemma 2.1) is un- 
affected by multiplication by holomorphic monomials. 

Hence, ak = Pn,mAk + Pn,mBk is orthogonal to the intersection of 
these sets; namely. 



{z^w^ : ^ {n,k) and {i^j) ^ (n,m)}. 



(4.1) 



Since 



^k-i -L {z^w^ : ^ {n,m — k — 1) and (i, j) ^ (n, m)} 



and since ak = am-k-i = z'^^'w'"' ^am_fc_i(l/z, 1/w), 

ak : ^ {n,m-k-l) and (z,j) ^ (n,m)} 

= {^V : {n,k) i and (n, -1) ^ (^,j)}. (4.2) 



Hence, is orthogonal to the union of the sets in (4.1) and (4.2). 
The set in (4.2) contains {z^w^ '■ i < n,j > 0} and the set in (4.1) 
contains {z'^w^ : i > n, j < m — 1}. Also, the set in (4.1) contains 
{z"'w^ : k < j < m — 1} while the set in (4.2) contains {z'^w^ : < j < 
k}. Combining all of this we get ak -L Ok- 

Finally, L is orthogonal to the intersection of Oq, . . . , Om-i- D 

We now look at the space generated by shifting the a^'s by powers 
of z. 

Theorem 4.1. With respect to L^( , '^'^,2 ), 

^ \Pn,m. I ' 

span{z'^aj{z,w) : < i, < j < m} 

= span{z^w-' : < z, < j < m} span{z^w-' : < i < n,0 < j < m} 

(4.3) 

and this is orthogonal to the larger set 

spa.n{z'^w^ : i < n, j > 0}. 
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Proof. Since the Ofc are polynomials of degree at most m — 1 in w, it is 
clear that 

spaji{z^aj{z, w) : < i,0 < j < m} C spanj^^w-^ : < < j < m}. 



By Theorem |3.1[ the are orthogonal to the spaces 

spaTi{z'^w^ : i < n, j > 0} D spa.n{z^w-' : i < n,0 < j < m}, 

and since these spaces are invariant under multiplication by z, the 
polynomials z^ak are also orthogonal to these spaces for all i > 0. So, 

spa,n{z^aj{z, w) : < i,0 < j < m} ± spa.n{z^w^ '■ i < n, j > 0}. 

Therefore, 

span{z'' aj{z, w) < k,0 < j < m} 

C spa.n{z^w^ '■ < i,0 < j < m} spa.ia{z^w^ : < i < n,0 < j < m} 

(4.4) 

and this containment must in fact be an equality. 
Indeed, any / in 

span{2;*w-' : < i,0 < j < m} 

which is orthogonal to {z''aj{z,w) : < k,0 < j < m} satisfies 
f{i,j) = for i > 72 and < j < m by Remark 3.3 and equation 
(3.2). Such an / cannot also be orthogonal to the space spaii{z^w^ : 
< i < n,0 < j < m} without being identically zero. □ 

Define 

H =span{2;*w-' : (0,0) < < {n,m- 1)} 

e span{;zV : (0,0) < (i,j) < {n-l,m- 1)}. 

Define also the reflection H 

H =span{zW : (0,0) < (i,j) <{n,m- 1)} 

Q spa.n{z''w^ : (1,0) < (z,j) < {n,m- 1)}. 

Proposition 4.2. We have the following orthogonal direct sum decom- 
positions in L^{l/\pn,m\'^da) 

oo 

■.0<k,0<j<m} = ^z'H (4.5) 

1=0 

oo 

"Hi := spMi{z^w^ : < A;, < j < m} = z'H. (4.6) 



2=0 
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If Kh is the reproducing kernel for H and is the reproducing kernel 



for H, then the reproducing kernel for the spaces in (4.5) and (4.6) are 

given by 



ZZi 



respectively. 

Proof. Now if is an m dimensional space of polynomials contained in 
the space (4.3) of the previous theorem. In particular, 

± span{zW : 2 < n,j > 0}, (4.7) 

and 

H = span{^*w-^ : ^ < < J < m} spa.n{z'^w^ : i < n,0 < j < m} 

since this space is also m dimensional and contains H. From this it is 
clear that H X z^H for z > and we have 



^ z^^H =spa.n{z^w^ : < j < m} 



i=0 



Q span{z'^w^ : i < n,0 < j < m}. 



Since shifts of H are contained in span{z'^w^ : < i, < j < m}, we 
must have 

oo 

^ z'H =spaii{z'w^ : < i, < j < m} 

1=0 

Q spaji{z'^w^ : 0<i<n,0<j< m} 

which combined with (4.3) gives (4.5). 

Next, H is also m dimensional and by (4.7) is orthogonal to 

{z^w^ : i > 0;j < m} 

which in particular contains the strip {z'^w^ : i > 0;0 < j < m} = 
z{z^w^ ■ i > 0;0 < j < m}. So, 

H =span{z^w-' '■ < i,0 < j < m} 

Q z spsiii{z^w^ '■ < i,0 < j < m} 
by dimensional considerations. Therefore, 

n^ = ^z^H. 

The formulas for the reproducing kernels are direct consequences of 
the orthogonal decompositions (see |1] for more on this). □ 
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Lemma 4.3. In L'^(l/\pn,rn\^da) the reproducing kernel for 
U = spmj^V : (0,0) < {i,j) ^ {n,m)} 

is 



Pn,m{z,w)pn,m{zi,Wi) - l)^{z, w)l)^{zi, lUi) 



(1 — ZZi){l — WWi) 



Proof. First, 



K{z,w;zi,wi) = K(^,^^^^){z,w) 



Pn 



[z,w)pn,m{Zl,Wi) 



(1 — ZZi){l — WWi) 

is the reproducing kernel for span{z'^w^ : (0,0) < {i,j)} since 



f{z,w) 



12 Pn,m[Z,W 
fizi,Wi 



:Pn,m{Zi,Wi)- 



dzdw 



{2T[iYz'w{l — zzi){l — wwi) 

Pn,m{zi,Wi) = f{zi,Wi) 



Pn,ni[Zl,Wi) 

by the Cauchy integral formula. On the other hand, 

Pn,m{z,w)^:^{Zi,Wi) 
(1 — ZZi){l — WWi) 

is the reproducing kernel for 

spm{zV : (0,0) < {i,j)}en. 

To see this it is enough to show that {z^w^pn,m '■ ( 0, 0) < {i,j)} is 
an orthonormal basis for the space (4.9). By Lemma 2.1, z'^w^pn,m is 



(4.8) 



(4.9) 



in the space in (4.9) for every z,j > and it is easy to check that 



these polynomials form an orthonormal set. We show that their span 
is dense. 

We may write Pn,m = cz"'w"^ + lower order terms with c 7^ 0, since 
Pn,m is stable. Now, let / be in the space in (4.9). If / ± Pn^m = 
cz'^w"^ + lower order terms, then since / is already orthogonal to the 
"lower order terms" we see that / ± z^w"^. Inductively, then, we see 



for alli< N and j < M but (z, j) ^ {N, M) 



that assuming / 
and assuming / 

/ will be orthogonal to the lower order terms in z^ w'^^pn^m- Therefore 



z W-' 

z^w^'^Pn,m, we automatically get / 



if / in (4.9) is orthogonal to {z^w^pn,m '■ i,j > 0} there can be no 
minimal {i,j) > {n,m) (in the partial order on pairs) such that / is 
not orthogonal to z'^wK In particular, / ± z'^w^ for alH > n and j > m 
and by ( 4.9[ ) / ± "H, which forces / = 0. 

So, {z'^w^Pn^m '■ (0,0) < (i,j)} is an orthonormal basis for the space 
in ( |4.9 ) while (4.8) is the reproducing kernel for this space. 
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Finally, the reproducing kernel for 

H = span{z'w^ : (0,0) < e (span{;2V : (0,0) < ihj)}en) 

is the difference of the reproducing kernels we have just calculated. 
Namely, 

Pn,m{z, W)pn,m{zi, Wi) - Ih^jz, w)l)^{zi,~Wi) 
(1 — ZZi){l — WWi) 

□ 

5. The bivariate Christoffel-Darboux formula 

Set 

Hi =spajLi{z^w^ :0<i<n,0<j<m — 1} 

spa.n{z^w-' : < i < n — 1,0 < j < m — 1} 

and 

H2 =span{2;*w-^ :0<i<n — l,0<j< m} 
span{z^w^ , < i < n — 1,1 < j < m}. 
The two variable Christoffel-Darboux formula is the following. 

Theorem 5.1. Let pn.m be a stable polynomial. Let Ki he the repro- 
ducing kernel for Hi and let K2 be the reproducing kernel for H2. Then 

Pn,m{z, W)pn,mizi, Wi) - ^^{z, w)^Z^{Zi,Wi) 

= (1 - wivi)Ki{z,w; zi,wi) + (1 - zzi)K2{z,w; zi,wi). 
Proof. Set 

n = spCnj^V : (0,0) < {i,j) ^ {n,m)} 
< i,0 < j < m} 
<i <n,0 < j} 



Hi = spaTi{z'^w^ 

= spanl^^w-^ 

and notice that "Hi and together span "H. 
Theorem says 

span{z*aj(z, w) : < i,0 < j < m} (5.1) 

= nie{'Hin'H2) cnen2 (5.2) 

which a fortiori implies 

'Hie(ninn2) = nen2. 

To see this, suppose / G {HeH2)e{Hie{ninH2)). Then, / e HeHi. 
As Til and 7^2 span Ti, such an / must be orthogonal to all of Ti and 
must equal 0. 
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The reproducing kernel for the space 

Ki(z,w;zi,wi) 



H e = z^Hi is 



1 — ZZi 



from Proposition 4.2 



If we interchange the roles of z and w in Proposition 4^ we see that 

K2{z,w; zi,wi) 
1 — WWi 

is the reproducing kernel for 'H2- 

Finally, the reproducing kernel for "H = ©("H 0^/2) can be written 
in two ways. On the one hand it equals 

Pn,m{z, W)pn,m{zi, Wi) - thZi^jz, 

(1 — ZZi){\ — WWi) ' 

but on the other it equals 

K2{Z,W] Zi.Wi) ^ Ki{z,W]Zi,Wi) 
1 — WWi 1 — ZZi 

by the discussion above. Equating these formulas and multiplying 
through by (1 — 2:^1) (1 — wwi), yields the desired formula. □ 

6. Parametric orthogonal polynomials 

The above results also shed light on the parametric orthogonal poly- 
nomials. The following proposition shows that the inner products of 
ao, . . . , ttm-i with respect to L^{dfi^, T) for the measures parametrized 

hj z = e'^ eT 

^ ^ ' 27r|p„,^(e^^w;)|2 ^ ' 

are trigonometric polynomials in z. 

Proposition 6.1. For fixed z eT 

• Mz,w;v)\ \ '^7' =-z-L{z,v;v) (6.2) 



and as a consequence 



\dw\ 



a.(z,t.)a,(.,z.)2^l^^^^^^^^^|,=T,,(.). (6.3) 



Proof. For 2; G T the expression 



-nrf N _ Pn,miz, w)pn,m{.Z, T]) - Pn,m{.Z, w)pn,m{.Z, T]) 

z L[z,w,ri) — z 

1 — wrj 
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is the reproducing kernel/Christoffel-Darboux kernel for polynomials in 
w of degree at most m—1 with respect to the measure \dw\/ {27i\pn,m(yZ, 
Indeed, this is one of the main consequences of the Christoffel-Darboux 
formula in one variable (see [7| equation (34) or [H] Theorem 2.2.7). It 
is a general fact about reproducing kernels K{w,ri) = Kr^iw) that 

\\K{;r^)\\^ = {K,,K,) = K{^,r^). 



Using these two observations, (6.2) follows. Equation (6.3) follows from 



matching the coefficients of rj'^rj^ in (6.2). □ 



Given Tm{z) defined in equation (2.1), set Di{6) as the determinant 



of the ixi submatrix of Tm{e ) obtained by keeping the first i rows and 
columns and set Do = 1 . We now perform the LU decomposition of 
Tjn which because it is positive definite does not require any pivoting. 
Set 

[<PL-iH, . . . , ^l{w)f = U{9)[w-\ . . . , 1]^, (6.4) 

where U{d) is the upper triangular factor obtained from the LU de- 
composition of Tm without pivoting. We find: 

Proposition 6.2. Suppose pn,m is a stable polynomial then {(j)f{w)}^^ 
satisfy the relations 

• 4>i{w) is a polynomial in w of degree i with leading coefficient, 
Dm-i+iiey 

which uniquely specify the polynomials. The above implies 



Proof. From the definition of we see that it is the inverse of the m x 
m moment matrix associated with dfi^{w). The first part of the result 
now follows from the one dimensional theory of polynomials orthogonal 
on the unit circle. The second part follows since z^^™''^^ Dm-j{0) is 
polynomial in 2;. □ 
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